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We study the polarization vector of t~ ( r + ) leptons created in v T (y T ) - Nucleon deep in- 
elastic scattering. Our work is oriented toward ultra-high energy particles relevant for cosmic-rays 
experiments. The goal of this paper is to derive relevant informations to constrain the systematic 
uncertainties that follow from the lack of knowledge about this point particularly when studying 
Earth-skimming neutrinos. 
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I. INTRODUCTION 

Through the last years, the observation of ultra-high 
energy (UHE) neutrinos has become one of the challenges 
of astroparticle physics. Many models, either astrophys- 
ical or exotic models, predict a substantial flux of neutri- 
nos. One of the most certain contribution to this neutrino 
flux are the so-called GZK-neutrinos produced in the de- 
cay of pions and kaons, from the interaction of ultra-high 
energy protons with the CMB [l[ . 

Such a mechanism provides a substantial flux of muon 
and electron neutrinos at the point of interaction. But 
given the large distances traveled by the particles, an 
observer can expect equal fluxes of v e , and v T at the 
observation point, due to flavour mixing and oscillations 
@,|. 

During the last years, an increasing effort has been 
put forward to develop a new generation of dedicated 
neutrino telescopes, both in the southern d, [j| and the 
northen @, 0, H, 0] hemisphere. Such detectors are rele- 
vant for an energy range of 1CU 6 EeV < E v < 1CU 1 EeV. 
At higher energies other experiments develop the detec- 
tion of coherent radio emission pro duced by neutrino- 
induced showers in matter [13, [Tllll2|| . Ultra-high energy 
cosmic-ray detectors such as the Pierre Auger Observa- 
tory [TU and the HiRes Fly's Eye detector [lj] , although 
they were not developed for the detection of neutrinos, 
may have equal or even better potential in the ultra-high 
energy range of 10 _1 EeV < E v < 10 2 EeV, where the 
GZK-neutrinos are expected 

In addition to the classical way of detecting neutrinos, 
at large zenith angles (8 > 75°) [H, it has been 
pointed out recently that the detection potential could 
be enhanced by the presence of v T , due to oscillations, 
in the cosmic neutrinos flux 

HIIIU. Upward- going UHE 
tau neutrinos that graze the Earth just below the horizon 
(also called "Earth-skimming neutrinos") have a quite 
high probability to interact in the crust and produce a tau 
lcpton which, if produced close enough from the surface, 
may emerge and trigger an extensive air shower which 
may be detected by a surface detector, provided it does 
not decay too far from the ground. 

The estimation of the sensitivity to such UHE neutri- 
nos implies the use of Monte Carlo simulations in which 
the propagation of tau neutrinos and tau leptons, pro- 



duced in v T charged current (CC) interactions, is simu- 
lated. To obtain the best description, all relevant pro- 
cesses for the energy scale considered must be taken into 
account. The ultra- high energy cosmic ray range dealt 
with here implies large extrapolations for the physical 
quantities considered in the calculation such as the par- 
ton distribution functions (pdf) involved in the simula- 
tion of the interactions with the nucleons of the crust. 
These large extrapolations lead to large systematic un- 
certainties which must be studied precisely to estimate 
the relevance of the different results obtained from these 
Monte Carlo calculations. In the estimate of the sensitiv- 
ity to Earth-skimming neutrinos, the uncertainties follow 
for a big part from the CC and NC deep inelastic cross- 
sections (related to the pdf uncertainties) and from the 
tau lepton energy loss along its propagation. Further- 
more, during v r CC interaction, the tau leptons created 
are likely to be highly polarized. This polarization plays 
also an important role in estimating the systematics as it 
can impact, mainly on the decay of the r lepton, and in 
a smaller extent on its propagation. When a tau decays, 
the energy distribution of the different products depends 
drastically on the helicity of the decaying tau lepton [2(| 
and this affects directly the development of extensive air 
showers (EAS) triggered by such decays, thereby influ- 
encing the detection and identification of such air show- 
ers. For surface detectors based on the study of EAS, the 
tau polarization can thus be one of the largest source of 
systematic uncertainties, through the calculation of the 
acceptance [2lj . 



The aim of this paper is to study more precisely the 
polarization of tau leptons produced in ultra-high energy 
v r CC interactions. Our goal is to constrain as much 
as possible this source of systematic uncertainties for ex- 
periments that are sensitive to the decay products of the 
tau lepton. We use a Monte Carlo calculation in which 
CC interactions of UHE tau neutrinos are simulated and 
the polarization of the induced tau lepton is derived for 
each interaction. In section [Til we present the different 
ingredients used to derive the tau polarization. In section 
III! we apply this framework to mono-energetic beams of 
tau neutrinos and present the results obtained. Finally, 
we discuss our conclusions in section ITVl 
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FIG. 1: Scheme of a nucleon deep inelastic interaction in 
the laboratory frame. 



II. CALCULATION OF TAU POLARIZATION 

We study the tau polarization through a Monte Carlo 
procedure. We simulate a neutrino beam that is forced to 
interact through CC interaction and for each interaction 
we compute the spin polarization of the produced tau. 

The reaction to be described is the following: 
v T (p v )N(p) — > t{p t )X, where q = p v — p T is the trans- 
fered 4-momentum. The x axis is defined along the di- 
rection of propagation of the r lepton. The y axis is 
chosen so that the scattering plane corresponds to the 
x — y plane. 

In all the following we consider an isoscalar nucleon 
target. 



A. Ingredients 

For the computation of the polarization, we followed 
the work of Levi [23| who derived practical formulae for 
the polarization vector of the tau lepton produced in deep 
inelastic neutrino-nucleon scattering. Once we have a 
simple formula to compute what we want to study (see 
equations [5J [9] and LTUJ) , we still need some other ingre- 
dients. First, simulating a CC deep-inelastic interaction 
is nothing more than drawing a pair of bjorken variables 
(x, y) that fully characterise our interaction. In the deep 
inelastic frame, these dimensionless variables are defined 
as follow 



E V ~E T 
E v 

Q 2 

2ME v y 



(1) 



These definitions are expressed here in the laboratory 
frame (fixed nucleon target) for simplicity but their val- 
ues are frame independent. E v is the incident neutrino 
energy E T the energy of the produced tau and Q 2 = —q 2 . 
M is the nucleon mass. We thus need an (x, y) distribu- 
tion to simulate the interaction. 

For the computation of these distributions as well as 
that of the tau polarization, we also need to use a pdf set 



to evaluate the distribution of partons inside the nucleon. 
Two different sets were used to check the dependency of 
our results upon this choice. We chose CTEQ6 [H| and 
MRST distributions [25| . that are both fitted to experi- 
mental data, where available. 

For the high energies considered in our work, regions 
of very small x can be probed during the scattering (see 
equation [5]) , where no experimental informations exist. 
We thus have to control the extrapolation to these small 
x in a secure way. For that we used the prescription of 
Reno [26] : 



xq(x,Q 2 ) = (-^— ") xq(x min ,Q 2 ), 



for x < x„ 



(2) 

We chose x m i n = 10 H , which is the lower limit for 
CTEQ6 set. A is determined for each flavour from the fit 
of the density function with Q = Myv . 

We derived our own distributions in x and y from the 
expression of the differential cross-section given in refer- 
ence 
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2G%ME V 



dxdy it 
x[xQ(x,Q 2 ) 



Ml 



M 2 W 



,Q 2 - 

xQ(x,Q 2 )(l-y 2 )}, 



(3) 

1.16632 x 



My/ is the mass of the W boson and Gf 
10~ 5 GeV -2 is Fermi's constant. 

Q and Q are the quarks distribution functions and are 
written from the individual distributions of valence and 
sea quarks obtained from the pdf sets mentionned above. 
For exact expressions see the previous reference. 

The formula for V T scattering is obtained from equa- 
tion [3] mostly by interchanging Q and Q, though it must 
be noted that the exact expressions for these two dis- 
tribution functions are not exactly the same for v T and 
v T . 

The physical regions for x and y are obtained by Al- 
bright and Jarlskog [28], with m T the mass of the t: 



< x < 1 



and 

where 

A = i ( i - 



2M(E U - m T ) 
A — B < y < A + B 



2ME„ 



2E l)i 



1 _l_ xM 

+ 2E V 



1 



2ME v x 
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(4) 
(5) 

(6) 



To the previous region for y given in equations [5] and 
O we added the condition: 



(Mj-M 2 ) 
2ME V 



<y<i 



E„ 



(7) 
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FIG. 2: Distribution in the log x — log y plane for the physical 
region of the v T — N charged-current interaction for an initial 
neutrino energy E„ = 10 GeV. We plot the entire available 
region. The pdf set used here is CTEQ6. The other set does 
not produce clearly visible differences. 



The upper limit is just a safety condition asking for the 
minimum energy of the produced tau to be m r (r pro- 
duced at rest). The minimum limit chosen here, with Ma 
the mass of the A resonance of the nucleon, requires that 
the interaction always occurs in the deep- inelastic region. 
This is convenient both to reduce the computation time 
and increasing the precision of our (x,y) distributions. 
For the high energies considered here, this is well justified 
as deep-inelastic scattering is totally dominant compared 
to quasi-elastic scattering or resonance production [29j . 
FigE] shows an example of distribution used for our cal- 
culation. It gives the physical region in the log(x)-log(y) 
plane, as defined above, for CC interaction of an incident 
v T of energy E v = 10 11 GeV. 

The dependance of such a distribution on the energy 
is quite easy to sum up. Increasing energy extends the 
physical region accessible in the x — y plane and shifts 
the pattern to smaller x. Of course, we use different 
distributions for v T and v r scattering. 



B. Tau polarization 

Once we have fulfilled all those requirements we can 
head to the calculation of tau polarization. The formula 
for the t~ lepton spin polarization quoted from the above 
reference reads: 



ME V 



(jt + ^)Q(x,Q 2 ) 

X 



7?(l 



--y)QM 2 ) 



(8) 



For the r + lepton produced in an anti-neutrino scat- 
tering we have: 



Sr 



ME,, 



(7? + ^)Q(x,Q 2 ) 
x 



--y)Q(x,Q 2 ) 

X 



Pq+P, 



(9) 



°- 2 In these expressions, Q and Q are defined exactly as 
above. S T = 2 me and 



Pr, 



X 



% = (1 -»)(!-!/ ' 



X 



(10) 



The 3- vectors ~jf and if are the momenta of the proton 
and of the transfered momentum q respectively. They 
are both expressed in the tau rest frame. To obtain the 
expression of the spin polarization in another frame, one 
just have to boost the 4-vector (0, ~P). 

In the laboratory frame, the spin polarization vector 
always lies in the scattering plane [29| so that the spin 
polarization vector in the tau rest frame can be written 
as 

~P = (P x ,P y ,P z ) = P(cos0p,sinflpcos0 P ,O), (11) 

where 8p and <f>p are the polar and azimuthal angles of 
the spin vector in the tau rest frame, respectively, and P 
denotes the degree of polarization. P = 1 corresponds to 
a fully polarized tau, whereas P — gives the unpolar- 
ized case. 

In our case, the useful physical information lies mostly 
in P x /P = cos Op. The massless particle case would lead 
to P x = ±1, which means fully left-handed t~ or right- 
handed t + . For the non-zero mass r leptons this is not 
the case anymore. 

In the following, we will be interested in evaluating 
the value of P x /P, more particularly in the laboratory 
frame for the UHE range of KT 1 EeV <E V < 10 2 EeV. 
This range is relevant for surface detectors such as the 
Pierre Auger Observatory and corresponds to the detec- 
tion window for the expected GZK neutrinos. For such 
experiments, the tau polarization represents a large sys- 
tematic error. This is easy to understand. Depending on 
the created t^ 1 helicity, the decay products have differ- 
ent energy distributions. This influences the efficiency of 
trigger and of identification and thus the acceptance of 
the detector to the emerging tau leptons [U [22j . From 
the previous references, we quote a 30% difference for 
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the acceptance to between the two extreme cases of 
polarization, i.e. all —1 or all +1 helicities. It must be 
noted that r + with +1 helicity gives the same energy dis- 
tributions than —1 t~ , so that they should lead to the 
same acceptance. This symmetry is nearly exact, even if 
small corrections remain in the angular distributions of 
the multihadron final states [2(3], sensitive to CP parity. 
But this effect is of course negligible and a full symmetry 
can be assumed between r~ and t + . 

Our study aims at improving our knowledge on the 
helicity of the produced leptons in order to constrain 
the calculation of the acceptance to these particles and 
thus to reduce the systematic error due to the lack of 
information on the polarization. 

We will also derive some informations for lower ener- 
gies, which could be useful for experiments working in 
this lower range of energy. 



III. RESULTS 

We simulated mono-energetic beams of v T and V T 
which interact through CC interaction with an isoscalar 
nucleon at rest, in other words, cosmic (anti-)neutrinos 
that interact with the Earth crust. 

We are interested in studying different features of the 
tau polarization both in the tau rest frame and the lab- 
oratory frame: the evolution with energy, differences be- 
tween t~ and r + leptons and the influence of different 
pdf sets on our results. 

For the discrepancy due to the different pdfs, if we 
expect a difference, it is at high energy. In figure [3] we 
show the polarization degree P for leptons induced 
by incident neutrinos of 10 4 and 10 11 GeV. This value 
is given here in the tau rest frame. In red, we plot the 
distribution obtained using CTEQ6 set and in blue that 
from MRST. 

The result obtained is in accordance with what we ex- 
pect: no differences at low energies, contrary to what we 
observe at higher energies, where existing informations 
constraining the values of pdfs have to be extrapolated 
to regions where we have no way to evaluate the physical 
evolution. An other information is the evolution of the 
polarization degree of the produced t~ leptons. The cre- 
ated particles are highly polarized even at low energies. 
At first order, one could almost consider the fully polar- 
ized case, whatever the energy. For produced r + leptons, 
the same comments hold. 

In figure [U we now show the distributions obtained for 
P x /P, namely the helicity, in the laboratory frame, for 
E v = 10 4 GeV to 10 7 GeV. The red histograms are for 
t~ and the blue ones for t + . For details on the figures, 
see the caption. The differences between the tau leptons 
and their antiparticles appear clearly. They should be 
reduced as energy increases because the parton distribu- 
tions involved in both interaction (v T N and V T N) be- 
come more and more dominated by sea quarks and tend 
to the same value. We do not talk here about the differ- 



| Mean 0.9995 ± 2.786e-oF 




0.75 0.8 0.85 0.9 0.95 1 

Polarization degree 

Mean 0.9999 ± 8.68e-07~ 
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FIG. 3: Degree of polarization for r~ leptons produced in v T 
CC interaction, from a mono-energetic v T beam. The upper 
figure corresponds to initial neutrinos of energy E v = 10 4 
GeV and the other to E v = 10 11 GeV 



ence in the overall sign, which is obvious, but more about 
the difference in the shape of the distributions which in 
consequence should get almost symmetric at very high 
energy. The r + leptons also tend to be slightly more 
polarized than the r~ at low energy. 

It is obvious that the produced leptons get more and 
more polarized as the energy increases. For r~ , the dis- 
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p x /p p x /p 




FIG. 4: Helicity of r leptons produced in CC interactions. 
In each graph, we plot in red the distribution for r~ leptons 
and in blue the one for r + leptons. In the upper figure, the 
solid line corresponds to the distribution for E v = 10 4 GeV 
and the dashed line to E v = 10 GeV. In the bottom figure, 
solid line is for E u = 10 6 GeV and dashed line for E v = 10 7 
GeV 



tribution seems to tend to only two bins, —1 and +1, the 
proportion of —1 helicity leptons growing with energy. 
For t + the same conclusions hold by inverting —1 and 
+1 helicities. 



FIG. 5: Helicity of r leptons produced in CC interactions. 
We plot in red the distribution for r~ and in blue the one for 
t . We consider incident neutrinos of energy E v = 10 8 GeV 



At higher energies, in the region that we are interested 
in, this behaviour becomes still more obvious. We show 
in figure [5] the same distribution for incident neutrinos of 
energy E v — 10 s GeV. We use the same color code than 
above. 

The distribution consists now only in 2 bins in —1 and 
+1, which means only left-handed or right-handed parti- 
cles. We also computed the distributions for higher ener- 
gies but it is not useful to show them here as they do not 
give us much more information than what we learn from 
the previous plot. All that we see is that every created r~ 
(resp. r + ) are mostly left-handed (resp. right-handed) 
with a ratio of 1000 : 1 in favour of —1 helicity (resp. 
+1). At E v = 10 11 GeV the mean helicity is -0.9983 
(resp. 0.9982). 

This result is exactly what is needed in order to address 
the problem of the systematic uncertainties due to the 
tau polarization. In view of the values shown above, it 
is totally fair to consider that all produced r~ have — 1 
helicity and that all r + are produced with a +1 helicity. 

All the previous results were obtained by using CTEQ6 
pdf set. We also computed the same distributions using 
MRST to check the influence of the set used. At high- 
est energies, we obtain a slightly larger proportion of +1 
and —1 helicity for r~ and t + respectively but this is 
only a second order effect. The two different sets used 
in this work represent well the range of the distribution 
functions present throughout the litterature and hence 
the present result can reasonably be assumed to be inde- 
pendant from the pdf choice. 
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IV. CONCLUSION 

In this paper we have studied the polarization vector of 
leptons produced in v T or V T deep- inelastic interaction 
with the Earth crust, in order to address the problem of 
the systematic uncertainties in the frame of large surface 
detectors such as the Pierre Auger Observatory. 

Our main result allows to choose unambiguously be- 
tween — 1 or +1 hclicities, for both r~ and r + leptons, 
when studying their propagation or decay. This result 
has been shown to be only lightly dependant on the par- 
ton distribution function set used. The ratio — 1/ + 1 
particles may be slightly affected, but this is just a sec- 
ond order effect that cannot impact deeply on any study 
following our conclusions. 

The last point, along with the fact that —1 helicity 
and +1 helicity r + lead to the same energy distributions 
among their decay products, allows to choose between 
the two extreme cases of acceptance to up-coming 
leptons and so to eliminate the r polarization as a source 



of systematic uncertainties. 

This conclusion holds when considering usual pdf sets 
but the study could be repeated using exotic models for 
the evolution of the pdfs, to quantify the influence of new 
physics. 

Below 10 8 GeV, even if the —1 helicity case is prevail- 
ing, the picture is not so contrasted, so that for studies 
concerning this lower energy range, extra care should be 
taken in evaluating the tau polarization distribution and 
its impact on their particular problem. 
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